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ABSTRACT 

Degradation  products  that  could  reasonably 
form  upon  storage  of  sulfur  mustard  have  been 
synthesized  and  analyzed  by  GC/MS. 
Elucidation  of  the  fragmentation  pathways  was 
accomplished  by  means  of  B3LYP/6-31+G(d) 
calculations. 

INTRODUCTION 

Sulfur  mustard,  HD  or  bis(2- 
chloroethyl)sulfide  1  was  first  introduced  as  a 
chemical  warfare  agent  in  WW 1 .  Its  production 
has  continued  to  the  present,  most  horrifically 
against  civilians  in  the  1980’s  Iran-Iraq  war. 
Due  to  its  ease  of  synthesis  and  stockpiling  by 
several  countries,  its  exploitation  by  terrorists 
remains  an  alarming  possibility.  Numerous 
publications  exist  concerning  the  pathways  and 
products  of  sulfur  mustard  degradation  under 
both  field  and  laboratory  conditions.  One 
compound,  2,  has  been  identified 
spectroscopically  in  ton  containers  in  quantities 
approaching  3%;  however,  direct  confirmation 


1  2 
of  its  presence  has  not  been  determined  using  an 
authentic  synthetic  sample.  In  order  to 
thoroughly  study  the  degradation  products  of  1, 
2  and  ten  related  sulfides,  sulfoxides  and 
sulfones  that  could  reasonably  form  in  the 
degradation  and  manufacture  of  1  were  prepared. 
These  compounds  had  the  following  structures: 


2-4,  R  =  CH2CH2C1 ,  x  =  0  -  2,  respectively; 
5-7,  R  =  CH2CH2OH,  x  =  0  -  2,  respectively; 
8  -  10,  R=  CH=CH2,  x  —  0  —  2,  respectively; 


11  and  12,  x  =  0,  and  R  =  CH2CH2OCH3  and 
CH2CH2OCH2CH3,  respectively.  Of  the  arsenal 
of  analytical  methods  available  for 
characterization  of  these  organosulfur 
compounds,  mass  spectrometry  (MS)  offered 
the  significant  advantages  of  establishing  the 
molecular  weights  of  the  compounds  (from  the 
molecular  ion  produced  by  abstraction  of  one 
electron),  and  in  addition,  generation  of 
fragmentation  patterns  that  are  exquisitely 
sensitive  to  specific  molecular  structure. 

While  it  is  possible  to  intuit  fragmentation 
patterns  from  available  empirical  MS  data,  it  is 
preferable  to  justify  mass  spectral  assignment  in 
terms  of  sound  theoretical  principles. 
Accordingly,  calculations  were  performed  on 
each  mass  spectral  peak.  Transition  states  for 
transformation  of  the  fragments  were  identified 
in  order  to  estimate  the  relative  intensity  of  each 
peak. 

Gas  chromatography-mass  spectrometry  was 
performed  on  these  organosulfur  compounds. 
Analysis  was  effected  both  with  electron 
ionization  (El)  and  collision-induced  dissociation 
(CID)  using  argon  as  the  neutral  gas.  The  El 
mass  spectrum  typically  contains  the  molecular 
ion  and  many  fragment  ions  providing  useful 
structural  information.  CID  is  a  technique  used 
to  confirm  the  structure  of  the  selected  mass  ion. 

The  synthesis  of  2  involved  the  nucleophilic 
attack  of  2-thioethanol  on  1  in  the  presence  of 
NaOCH3,  followed  by  conversion  of  the  terminal 
alcohol  functionality  into  a  chloro  substituent  in 
77%  overall  yield.  The  remaining  potential 
degradation  or  side  products  in  the  manufacture 
of  1  involved  oxidation  chemistry  with  H202, 
elimination  of  HC1  using  triethylamine,  or  KOH 
in  methanol  or  ethanol. 

The  nature  of  the  mass  spectral  fragments 
demands  calculations  at  a  reasonably  high  level 
that  incorporate  electron  correlation. 
Considering  the  computational  cost  and 
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anticipated  accuracy  of  the  results,  density 
functional  theoretical  (DFT)  methodology  was 
chosen.  (Bickelhaupt,  2000;  Cramer,  2004)  It 
was  decided  to  use  the  B3LYP/6-31+G(d)  level 
for  which  there  is  considerable  literature 
precedent  for  calculations  in  mass  spectral 
elucidation.  (Chu,  2008;  Baciocchi,  2008; 
Anand,  2004,  Pepe,  1999.) 

In  summary,  eleven  compounds  relevant  to  the 
degradation  of  1  upon  storage  were  synthesized 
and  analyzed  by  gas  chromatography/mass 
spectrometry.  The  fragmentation  patterns  were 
elucidated  by  means  of  high  level  DFT 
calculations. 


COMPUTATIONAL  METHODS 

Initial  geometries  for  the  ground  and 
transition  states  for  the  cations,  cation  radicals 
and  neutral  species  were  obtained  by  means  of 
semi-empirical  calculations  using  AMI 
methodology.  The  final  DFT  calculations  were 
performed  by  means  of  the  Gaussian  suite  of 
programs  (Frisch  et  al.,  2004)  at  the  B3LYP/ 
6-31+G(d)  level.  A  scaling  factor  (Foresman, 
1996)  of  0.9804  was  used  for  the  thermal 
correction  to  the  computed  energies  at  a 
temperature  of  298.1  K.  Criteria  for  the 
transition  and  ground  states  were  calculation  of 
one  and  zero  imaginary  frequencies, 
respectively.  Gausview  3.09  was  used  to 
animate  the  sole  imaginary  frequency  and 
confirm  that  the  transition  state  was  found  on  the 
potential  energy  surface  connecting  reactants  and 
products;  IRC  calculations  were  used  to  confirm 
identification  of  the  transition  state.  Transition 
state  theory  was  employed  in  calculation  of  the 
free  energies  of  activation.  (Lowry,  1981.) 

GAS  CHROMATOGRAPHY-MASS 
SPECTROMETRY 

An  Agilent  5975  mass  spectrometer  interfaced 
to  an  Agilent  6890  gas  chromatograph  was  used 
for  gas  chromatographic-mass  spectral  analysis. 
An  Agilent  J&W  Scientific  HP-5ms  bonded 
phase  capillary  column  (30  m  x  0.25  mm  i.d.) 
with  a  film  thickness  of  0.25  pm.  The  injection 
port  temperature  was  220  °C,  the  GC-MS  interface 
temperature  was  250  °C  and  the  source  temperature 
was  150  °C.  The  carrier  gas  was  helium  with  a  flow 
rate  of  1  mL/min  and  the  oven  temperature  was 
programmed  from  60  to  250  °C  at  15  °C/min.  A  split 
injector  was  used  (split  ratio  75:1)  and  0.2  pL  volume 


of  sample  was  placed  on  the  column.  The  mass  range 
scanned  was  from  50  to  450  Da  at  4  scans  per  second. 

The  gas  chromatography-tandem  mass 
spectrometry  (GC-MS-MS)  analysis  was  performed 
on  a  Finnigan  TSQ-7000  GC-MS-MS  (San  Jose,  CA, 
U.S.)  interfaced  to  a  Varian  3400  gas  chromatograph 
(Palo  Alto,  CA,  U.S.).  The  instrument  was  equipped 
with  a  30  m  x  0.25  mm  i.d.  HP-5ms  capillary  column 
with  a  film  thickness  of  0.25  pm.  The  injector 
port  temperature  was  250  °C,  the  GC-MS 
interface  temperature  was  250  °C  and  the  source 
temperature  was  150  °C.  The  carrier  gas  was 
helium  with  a  flow  rate  of  1  mL/min  and  the 
oven  temperature  was  programmed  from  60  to 
250  °C  at  15  °C/min.  A  split  injector  was  used 
(split  ratio  75:1)  and  0.5  pL  volume  of  sample 
was  placed  on  the  column  manually.  Argon  was 
used  as  the  collision  gas  for  the  CID.  The 
pressure  in  the  collision  cell  was  1.5  mT  and  the 
collision  energy  was  -15  eV.  The  mass  range 
scanned  was  from  20  to  200  Da  at  1  scan/s. 
Qualitative  and  quantitative  analyses  were 
processed  by  Xcalibur™  (Finnigan  MAT,  San 
Jose,  CA,  U.S.)  software  supplied  with  the  MS 
data  system.  The  mass  ion  of  interest  was 
selected  for  each  individual  compound  and  the 
retention  time,  peak  area  and  relative  percent 
abundance  of  the  fragment  ions  were  calculated 
(Table  1). 


SYNTHETIC  PROCEDURES 

The  compounds  studied  herein  were  prepared 
by  previously  reported  procedures  (Winemiller, 
2008.)  As  an  example,  the  preparation  of  2  was 
accomplished  by  the  nucleophilic  displacement 
of  chloride  from  bis(2-chloroethyl)ether  by  2- 
thioethanol  in  the  presence  of  NaOMe  followed 
by  chlorination  of  the  product  with  SOCl2. 

MASS  SPECTROMETRIC  RESULTS 

2:  The  El  mass  spectrum  exhibits  the  molecular 
ion  at  m/z  202  and  isotopic  ions  at  m/z  204  and 
m/z  206  due  to  chlorine.  Higher  mass  ions  due 
to  [M-HC1]+,  [M-0CH2CH2C1]+  and  [M- 

CH2SCH2CH2C1]+  are  observed  at  m/z  166,  m/z 
123  and  m/z  93.  The  mass  ion  at  m/z  123  was 
selected  for  CID  and  produced  three  major 
product  ions  at  m/z  63,  m/z  87  and  m/z  95  due  to 
[C2H4C1]+,  [C4H7S]+  and  [C2H4SC1]+ 

respectively,  with  a  relative  intensity  ratio  of 
100:1:13. 


3:  The  El  mass  spectrum  shows  three  principal 
ions  with  m/z  values  of  63,  76  and  107.  The 
higher  mass  ion  at  m/z  107  was  formed  due  to 
[M-C2H40SC1]+  and  the  lower  mass  ion  at  m/z 
63  was  from  loss  of  C4H8S02C1.  The  mass  ion  at 
m/z  107  was  selected  for  CID  and  produced  two 
major  product  ions  with  m/z  63,  [C1CH2CH2]+ 
and  m/z  45,  [C2H50]+  with  a  relative  intensity 
ratio  of  100:12  . 

4:  The  El  mass  spectrum  shows  five  principal 
ions  with  m/z  values  of  63,  106,  127,  155  and 
185.  The  higher  mass  ions  at  m/z  155  and  m/z 
185  were  formed  due  to  [M-C2H50C1]+  and  [M- 
CH2C1]+.  Two  higher  mass  El  fragmentation 
ions  (m/z  155  and  m/z  185)  were  selected  for 
CID.  Three  product  ions  at  m/z  63,  m/z  93  and 
m/z  127  due  to  [C1CH2CH2]+,  [C2H5S02]+  and 
[C2H4S02C1]+  respectively  with  a  relative 
intensity  ratio  of  100:8:24,  were  observed  in  the 
CID  spectrum  of  m/z  155.  The  CID  of  m/z  185 
contains  three  major  product  ions  with  m/z 
values  of  63  (loss  of  C3H6S03),  127  (loss  of 
C3H60)  and  155  (loss  of  CH20),  with  a  relative 
intensity  ratio  of  25:26: 100. 

5:  The  El  mass  spectrum  exhibits  a  molecular 
ion  at  m/z  184  and  an  isotopic  ion  at  m/z  186  due 
to  chlorine.  Higher  mass  ions  due  to  [M-HC1]+ , 
[M-0CH2CH2C1]+,  and  [M-0C3H6C1]+  are 
observed  at  m/z  148,  m/z  105  and  m/z  91 
respectively.  The  mass  ion  at  m/z  105  was 
selected  for  CID  and  produced  three  major 
product  ions  with  m/z  values  of  87,  61  and  45 
with  a  relative  intensity  ratio  of  13:42:100.  The 
m/z  87  was  due  to  loss  of  H20,  the  m/z  61  was 
due  to  loss  of  C2H40  and  the  m/z  45  was  due  to 
loss  of  CH2CH2S. 

6:  The  El  mass  spectrum  contains  the  molecular 
ion  at  m/z  200  and  an  isotopic  ion  at  m/z  202  due 
to  chlorine.  Higher  mass  ions  due  to  [M-HC1]+ , 
[M-C2H5S02]+  and  [M-0C5H9C1]+  are  observed 
at  m/z  165,  m/z  107  and  m/z  94  respectively.  The 
lower  mass  ion  at  m/z  63  was  due  to  loss  of 
C4H9S03.  The  mass  ion  at  m/z  107  was  selected 
for  CID  and  produced  two  major  product  ions  at 
m/z  63  and  m/z  45,  due  to  [C2H4C1]+  and 
[C2H50]+  respectively  with  a  relative  intensity 
ratio  of  100:12. 

7:  The  El  mass  spectrum  shows  six  principal 
ions  with  m/z  values  of  167,  137,  109,  93,  79  and 
63.  The  higher  mass  ions  at  m/z  167  and  m/z  109 
were  formed  due  to  [M-CH2C1]+  and  [M- 
CH2CH20CH2CH2C1]+.  Two  El  fragments  ( m/z 
137  and  m/z  167)  were  selected  for  CID.  Three 


product  ions  at  m/z  93,  m/z  109  and  m/z  119  due 
to  [C4H7S02]+,  [C2H5S03]+  and  [C2H5S02]+ 
respectively  with  a  relative  intensity  ratio  of 
8:40:1,  were  observed  in  the  CID  spectrum  of 
m/z  137.  The  CID  of  m/z  167  contains  three 
major  product  ions  with  m/z  values  of  45  (loss  of 
C3H6S03),  109  (loss  of  C3H60)  and  137  (loss  of 
HCHO),  with  a  relative  intensity  ratio  of 
15:34:100. 

8:  The  El  mass  spectrum  contains  the  molecular 
ion  at  m/z  166  and  isotopic  ions  at  m/z  168  and 
m/z  170  due  to  chlorine.  Higher  mass  ions  at 
m/z  131  and  m/z  106  were  formed  due  to  [M-C1]+ 
and  [M-C3H4S]+.  Three  potential  resonance 
structures  exist  from  the  loss  of  Cl  from  the 
molecular  ion.  The  lower  mass  ions  at  m/z  63 
(loss  of  C4H7SO)  and  m/z  43  (loss  of  CH2CH2C1) 
are  shown  in  Figure  8b.  The  CID  of  the  mass 
ion  at  m/z  166  depicts  three  product  ions  with 
m/z  values  of  43,  63  and  131  with  a  relative 
intensity  ratio  of  68:29:3. 


9:  The  El  mass  spectrum  of  8  contains  the 
molecular  ion  at  m/z  182  and  a  higher  mass  ion 
due  to  [M-C4H80C1]+  at  m/z  107.The  CID  of  m/z 
107  contains  two  major  product  ions  at  m/z  45 
and  m/z  63  due  to  [C2H50]+  and  [CH2CH2C1]+ 
respectively,  with  a  relative  intensity  ratio  of 
12:100. 

10:  The  El  mass  spectrum  shows  four  principal 
ions  with  m/z  values  of  63,  91,  106  and  149. 
Higher  mass  ions  at  m/z  149  and  m/z  106  were 
formed  due  to  [M-CH2C1]+  and  [M-C2H4S02]+. 
Lower  mass  ions  at  m/z  63  and  m/z  91  were 
formed  due  to  [M-C4H7S03]+  and  [M-C4H80C1]+ 
(Figure  10a).  The  CID  of  m/z  149  contains  two 
major  product  ions  at  m/z  91  and  m/z  119  due  to 
[C2H3S02]+  and  [C4H7S02]+  respectively,  with  a 
relative  intensity  ratio  of  100:99. 

11:  The  El  mass  spectrum  contains  the  molecular 
ion  at  m/z  198.  The  higher  mass  ion  at  m/z  119 
was  formed  due  to  [M-0CH2CH2C1]+.  Lower 
mass  ions  at  m/z  103  and  m/z  59  were  formed 
due  to  [M-0CH3-C1CH2]+  and  [M-C3H80- 
0CH2CH2C1]+.  The  CID  of  m/z  119  contains 
two  major  product  ions  at  m/z  59  and  m/z  87  due 
to  [C2H3S]+  and  [C4H7S]+  respectively,  with  a 
relative  intensity  ratio  of  100:1 1. 

12:  The  El  mass  spectrum  shows  molecular  ion 
at  m/z  212.  The  higher  mass  ion  at  m/z  111  was 
formed  due  to  [M-C1]+  forming  one  of  three 


resonance  structures.  Another  higher  mass  ion  at 
m/z  133  formed  due  to  [M  -  0CH2CH2C1]+  could 
be  from  one  of  two  resonance  structures.  The 
CID  of  m/z  133  contains  four  major  products  m/z 
45,  m/z  73,  m/z  87  and  m/z  105  due  to  [C2H50]+, 


RESULTS  AND  DISCUSSION 


Since  2  has  been  found  in  measurable 
quantities  in  mustard  storage  containers,  it  was 
decided  to  initially  focus  on  the  analysis  of  this 
compound  to  define  the  parameters  for  the 
remaining  compounds  in  the  series. 

The  El  mass  spectrum  of  2  contains  the 
molecular  ion  at  m/z  202  and  m/z  166,123  and 
93.  When  mass  ion  m/z  123  was  selected 
for  CID  three  major  product  ions  were  produced 
at  m/z  63,  m/z  87  and  m/z  95.  The  full  proposed 
fragmentation  pathway  is  illustrated  in  Scheme 
1. 


[C4H90]+,  [C4H7S]+,  and  [C4H9SO]+  respectively, 
with  a  relative  intensity  ratio  of  98:100:7:6. 
[C2H3S02]+  and  [C4H7S02]+  respectively,  with  a 
relative  intensity  ratio  of  100:99. 


In  order  to  verify  this  scheme,  DFT  calculations 
were  performed  on  the  molecular  ion,  and  also 
on  each  of  the  fragments.  In  addition,  the  free 
energies  of  activation  of  the  transition  states 
were  calculated  to  assess  the  plausibility  of  the 
intensity  of  the  fragment  ions.  The  results  of  this 
calculation  are  provided  in  Table  1.  In  order  to 
determine  the  appropriateness  of  the  level  of  the 
calculations,  a  higher  order  basis  set,  6-3 1 1+G(d) 
was  used  for  Reaction  C  and  the  result  compared 
to  that  obtained  using  our  standard  basis  set. 
Higher-order  basis  sets  demand  considerably 


SCHEME  1 


m/z  95 


m/z  87 


TABLE  1 


Reaction 

(see  chart) 

A 

B 

C 

D 

E 

F 

G 

C  [6-311++G(d)] 


[EcorriGS 

hartrees 

-1629.175406 

-1629.175406 

-1629.175406 

-1629.175406 

-1015.216590 

-1015.216590 

-1015.216590 

-1629.342507 


[EcorrllS 

hartrees 

-1629.175590 

-1629.094208 

-1629.118185 

-1629.113152 

-1015.098082 

-1015.160039 

-1015.126255 

-1629.282150 


AE*  AH* 
kcal/mol  kcal/mol 


-0.115 

0.478 

50.95 

51.55 

35.91 

36.50 

39.07 

39.66 

74.36 

74.96 

35.49 

36.08 

56.69 

57.28 

37.87 

38.47 

SGS 

SXs 

eu 

eu 

126.64 

117.20 

126.64 

124.73 

126.64 

127.85 

126.64 

119.16 

87.382 

90.713 

87.382 

88.052 

87.382 

87.184 

126.64 

127.85 

AS* 

AG* 

eu 

kcal/mol 

-9.433 

3.29 

-1.905 

52.1 

1.216 

36.1 

-7.476 

41.9 

3.331 

74.0 

0.67 

35.9 

-0.198 

57.3 

1.215 

38.1 

AG*  AG* 

kJ/mol  eV 

13.8  0.14 

218.0  2.26 

151.2  1.57 

175.3  1.82 

309.5  3.21 

150.1  1.56 

239.9  2.49 

159.4  1.65 


Higher-order  basis  increased  computational  time.  different  by  2.0  kcal/mol  providing  support  for 

The  free  energies  of  activation  were  only  our  computational  methodology. 


FIGURE  1.  TRANSITION  STATES  LEADING  TO  THE  FRAGMENTATION 
PATTERN  MOLECULAR  ION  OF  2. 
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State-of-the  art  calculations  using  this 
methodology  are  unable  to  predict  intensities  of 
the  fragment  ions  in  the  mass  spectrum. 
Obtaining  plausible  energies  and  structures  for 
the  stable  ions  and  the  transition  states  leading  to 
them  provides  strong  support  for  the  Scheme  1 . 


Table  2  provides  the  energies  of  the  fragment 
ions,  their  entropies  and  illustrates  that  all 
calculated  frequencies  are  positive  by  giving  the 
value  of  the  lowest  calculated  vibrational 
frequency. 


TABLE  2 


Energy 

Therm.  Corr. 

Corr.  Energy 

Entropy 

Lowest 

hartrees 

hartrees 

hartrees 

eu 

Frequency 

mz  87 

-554.5110071 

0.1053210 

-554.4077504 

77.7 

51.6 

mz  95 

-936.6790049 

0.0621410 

-936.6180819 

71.4 

225.0 

mz  123 

-1015.3364636 

0.1222700 

-1015.2165901 

87.4 

58.3 

mz  166 

-1168.5456462 

0.1742400 

-1168.3748213 

113.7 

25.8 

MZ  93 

-653.0224172 

0.0949950 

-652.9292841 

79.1 

79.1 

MZ  63 

-538.4656706 

0.0592530 

-538.4075790 

63.0 

476.2 

MW202 

-1629.3619685 

0.1902920 

-1629.1754062 

126.6 

11.5 

HC1 

-460.7979997 

0.0090240 

-460.7891526 

44.6 

2925.1 

thiirane 

-476.7874231 

0.0586370 

-476.7299354 

62.3 

629.4 

CH2SCH2CH2C1 

-976.2670189 

0.0895060 

-976.1792672 

85.5 

56.4 

ch2=ch2 

-78.5932689 

0.0541440 

-78.5401861 

53.7 

837.1 

och2ch2ci 

-613.9696705 

0.0691350 

-613.9018905 

71.9 

130.9 

och2ch2sch2ch2ci 

-1090.7882904 

0.1248430 

-1090.6658943 

100.7 

33.6 

CONCLUSIONS 

DFT  calculations  using  the  6-3 1G  +(d)  basis 
set  have  been  applied  to  the  elucidation  of  the  El 
and  CID  spectra  of  2,  found  in  HD  storage 
containers.  A  computation  using  a  higher,  more 
expensive  basis  produced  only  minor  changes  in 
the  energy  calculations.  While  computational 
reproduction  of  spectral  intensities  is  beyond  the 
capability  of  the  method  employed,  the 
calculations  were  able  to  establish  the  relative 
energies  of  the  radical  and  the  cationic  and 
radical  products,  and  establish  that  the  reactants 
and  products  could  be  connected  by  means  of  a 
reasonable  transition  state.  Our  conclusions  for 
additional  members  of  this  series  of  compounds 
have  been  similar  to  those  for  2. 
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